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Abstract 

It is shown that the transverse momentum distributions of particles 
emerging from the decay of statistical clusters, distributed according 
to a power law in their transverse energy, closely resembles that follow¬ 
ing from the Tsallis non-extensive statistical model. The experimental 
data are well reproduced with the cluster temperature T ~ 160 MeV. 

1. It is now well-documented that the transverse momentum distribution 
of various particles at high energy and in a very broad range of the trans¬ 
verse momentum, is correctly described by the Tsallis distribution. This was 
observed in all high-energy experiments IB-01. as well as in the recent phe¬ 
nomenological analyses 0 -im This observation is usually interpreted in 
terms of the statistical model of particle production, employing the Tsallis 
non-extensive statistics mm- Such interpretation!!, although very attrac¬ 
tive, meets a serious difficulty, however: it is indeed not easy to explain why 
any statistical model can apply at very large transverse momenta, where the 
perturbative QCD phenomena are known to dominate. 

This problem was recently addressed in a series of papers mm where 
the ideas derived from perturbative QCD (as applied to hard interactions), 
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1 The fundamental relation between the Boltzmann and Tsallis statistical models is 
clearly explained in H3J. 
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accompanied by the parton cascade (responsible for the jet fragmentation) 
were used to explain this puzzling result. 

In this note, following the general idea suggested in m ns] (c.f. also 
PI), we apply it to the statistical model of particle production which is 
rather successful in describing data on particle multiplicities (For a review, 
see, e.g. HU EE]). We show that the Tsallis formula can be recovered, 
to a good accuracy, in the model where the observed particles are decay 
products of clusters [19, 2Q\ [2lj which (i) decay according to the standard 
Boltzmann statistics and (ii) the distribution of their Lorentz factors follows 
a power law, as suggested by (perturbative) QCD. This observation indicates 
that the experimental validity of the Tsallis formula may be interpreted as 
another confirmation of the standard statistical model rather than that of its 
non-extensive Tsallis versioijE 

Indeed, the intriguing ’’unreasonable” success of the statistical model in 
description of multi-particle production in various processes and at various 
energies suggests that the final stage of the process of hadronization is dom¬ 
inated by the hadrons in the state of statistical equilibrium. It is also clear 
that the equilibrium cannot be global, as the observed spectra are far from 
isotropic. These observations lead naturally to the idea [19] that the transi¬ 
tion from the early state of the process, dominated by interactions between 
the hadronic constituents, most likely proceeds through an intermediate stage 
of clusters emitting the final hadrons according to the rules of statistical 
physics. 

If one admits that this process of cluster formation and thermal decay 
is a universal feature of hadronization, one is led to the conclusion that 
also the high transverse momentum jets hadronize in the same way (c.f. 
[2D]). It follows that the characteristic features of clustering should leave 
their imprints even in the region of hard physics. In the present paper we 
show that this picture, when combined with the power law distribution of 
the (transverse) Lorentz factor of the cluster, leads to transverse momentum 
distribution of the decay products which is very close to that of Tsallis (and 
thus also close to experiment). 

In the next section the idea of the statistical cluster is formulated and 
the transverse momentum distribution of its decay products is derived. The 
relation to the Tsallis distribution is discussed in Section 3. Summary and 

2 Another approach aiming at the explanation of the power law tails within statistical 
model is discussed in T2 1. 
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comments are given in the last section. 

2. Following the ideas explained above, the decay distribution of the 
statistical cluster at rest is taken in the form of the Boltzmann distribution 
which, for a cluster moving with the four-velocity u p becomes 

p(p\ u)d 2 p±dy = e~P Pf,u '‘d 2 p^dy ( 1 ) 


where (3 = 1 /T. 

Consider a cluster at rapidity Y moving in the transverse direction with 
the velocity v±. We have 

u 0 = \Jl + u\ cosh Y\ u z = sinh Y ; v z = tanli Y ; 

u± = 7Uj_; 7 = (1 — v 2 Y l/2 —>• \]l + 7 2 rF = 7. ( 2 ) 

The distribution of particle momentum is then 

p(p, y )dy = dyd 2 p J d 2 v ± dYG(v ± , Y)e~ p ^ m± cosh ( y -Y)-pp ± u ± c OS <p ( 3 ) 

where 0 is the angle between v and p± and where we have denoted 

7_l = yfl + u i = 7\/l — v l u x = sji[ ~ 1 (4) 

Integration over 0 and y gives the distribution of the transverse momentum: 

p(p)d 2 p = d 2 p j d 2 v ± dYG(v ± ,Y)K 0 [/3m ± ‘y A _] I o[l 3 * P± u ±\ ( 5 ) 

3. To evaluate the distribution of transverse momenta of the cluster decay 
products, one needs the distribution of the cluster transverse velocity v±. In 
this paper we study a power law in the transverse Lorentz factor 7 j_ (for a 
fixed cluster mass, this would correspond to a power law in its transverse 
energy). Thus we take 

d?v±dYG(v ± , Y) ~ G(Y)dY^l K d^ 1 _ (6) 

Given simplicity of this assumption, it was rather surprising to find that 
it leads to the distribution which closely resembles that of Tsalli.^, from 

3 Qualitatively, the result of this kind may be actually expected, as it is the well-known 

EH C2E1 [ 26 ] that the Tsallis formula is naturally obtained by adequate fluctuations of the 

parameters of the Boltzmann spectrum. 
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Figure 1: Transverse momentum distribution of pions, kaons and protons 
from the statistical cluster decay (dashed lines), normalized to 1 at p± = 0, 
compared to two Tsallis distributions (Eq.([7])) (full lines). 1 GeV < p± < 200 
GeV. T = 155 MeV, k= 6.5. Best fit from p_j_ = 0 to p± = 50 GeV. 

p± ps 100 MeV up to p±= 200 GeV. This was verified numerically for the 
cluster temperature in the region from 100 till 180 MeV and the power k from 
4 till 7, i.e. in the range covering the physical conditions one may expect in 
high-energy collisions.. 

An example of such calculation is shown in Figs [l]and [2] where the distri¬ 
butions of pions, kaons and protons evaluated using ([5]) and (ED with k = 6.5 
and T=155 MeV, are compared with the two versions of the Tsallis distri¬ 
bution [6j |8, 13], [23]: 

D\ = cmj_[l + (q- 1 )m_ L /T ts ] q/{1 ~ q) ; 

D 2 = c[l + {q-l)mjT ts ] 1 ^- ri \ (7) 
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Figure 2: Same as Fig{T]but for 0 < p± < 5 GeV. Lines: the Tsallis distri¬ 
bution. Crosses and stars: statistical clusters. 

where c is the normalization constant, q — 1 measures the deviation from the 
standard statistical model and T ts is the Tsallis temperature^]. 

One sees that, except at very small p±, below ~ 100 MeV, there is an ex¬ 
cellent agreement between the two formulations and for all kinds of particles. 
One also sees that for p± > 1 GeV it is difficult to distinguish between the 
two versions of the Tsallis distributions. For the distribution D { the Tsallis 
parameter T ts can be approximated by the simple relation T ts ~ (q — 1 )T. 
This is not true, however, for D 2 . In this case the relation between T ts and T 
is more complicated and, moreover, it depends substantially on the particle 
mass. 

Recently, a new analysis of transverse momentum distribution of charged 
particles in terms of the Tsallis distribution has been published [6j. To com- 

4 The form D\ is obtained by demanding maximum of the Tsallis entropy, i.e. thermo¬ 
dynamic equilibrium [24) 125]. The second form is the standard Tsallis distribution. 
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pare these results with our approach, we have evaluated the distribution 
following from the decay of a cluster for pions, kaons and protons and con¬ 
structed the distribution of charged particles, using the weights (1:1:2), as 
proposed in [6j. In Fig. [3] the results in the region from p _l = 0 till p± = 5 
GeV are compared with the Tsallis distribution from [6]. One sees that the 
agreement is very good, except at p± < 100 MeV. The parameters of the 
Tsallis distribution in this case are q — 1 — 0.150 and T ts —76 MeV, in good 
agreement with (6j. The region p± > 5 GeV is not shown because in this 
region one simply cannot distinguish between the two curves. 
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Figure 3: Transverse momentum distribution of charged particles from 
the statistical cluster decay (crosses), compared to the Tsallis distribution 
(dashed line) used in |6j (the first formula in (JTJ) ). 0 < p± < 5 GeV. T = 155 
MeV, k= 6.5. Best fit from p± = 0 to p± = 50 GeV. 

4. In summary, we have discussed the transverse momentum distributions 
of particles emitted in the decay of a statistical cluster. It was shown that if 


6 



the (transverse) Lorentz factor of the cluster follows a power law, the result¬ 
ing distribution is very close to that derived from the Tsallis non-extensive 
statistics. 

This result may be considered as a possible explanation of the surpris¬ 
ing observation that the Tsallis formula works not only at small transverse 
momenta (where the ideas of statistical equilibrium may be applicable) but 
even at transverse momenta as large as ~ 200 GeV. 

Some comments are in order. 

(i) It should be emphasized that the observed similarity between the Tsal¬ 
lis formula and that following from the statistical cluster decay, is only an 
approximation. Our results indicate, however, that it may be rather diffi¬ 
cult to distinguish experimentally between these two approaches. Perhaps 
the measurements at larger transverse momenta may be helpful, as the two 
distributions start to deviate from each other at energies above 200 GeV. 

(ii) We have been discussing emission of a single statistical cluster. As it 
is rather unlikely that a high-energy jet may fragment into a single cluster, 
production of many clusters must also be considered. Since our discussion 
concerns only the single-particle distribution, however, the results are insen¬ 
sitive to the number of clusters produced in a given event, provided they are 
emitted independently. 

(iii) Clearly, the power law assumed in ([6]) is only a phenomenological 
guess and should be treated as such. Its main advantage is the extreme 
simplicity (for more elaborate calculations see, e.g., [HJCGiJ .IB])- Needless to 
say, the parameter k remains free at the present stage, and cannot be reliably 
evaluated from theory. 

(iv) It has been shown recently [27] that the distributions of transverse 
momenta at various energies follow a scaling law, suggested by the saturation 
property of the parton distributions. An interpretation of this observation 
in terms of the Tsallis approach was proposed in [T6j [28] • It would be thus 
interesting to investigate how this scaling property of the spectra translates 
into the results shown in the present paper. 
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